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In measuring the ip(2S) radiative decays at BESII, contribution of the background is serious in 
most of the final states. To extract the number of signal events, a fit to the x 2 distribution of 
kinematic fit using signal and background components is proposed. Extensive Monte-Carlo simu- 
lations (MCS) are performed, and the results show that the shape of \ 2 distribution of the signal 
channel looks different from those of the background channels, thus it provides us by fitting the \ 2 
distribution of the data to extract the number of signal events. An input-output test is performed 
using MCS, and the uncertainty of the fit method is found to be less than 2%. 

PACS numbers: 02.50.Ng, 13.20.Gd, 29.90.+r 



I. INTRODUCTION 

The radiative decays of J/ip and if>', due to the di- 
rect coupling of outgoing photon to the charm quark, 
could be an excellent laboratory to study the hadron 
spectroscopy and search for exotic states. At the lead- 
ing order of pQCD picture, both the Jjif) and ip radia- 
tive decays can be described by the cc quark annihila- 
tion with the emission of one photon and two gluons. 
There could be a possibility that the two gluons might 
form a glueball, or one gluon combines with a qq pair 
to form a hybrid state, or the two gluons couple to 
a multiquark state. Compared with J/ij) decays, few 
measurements of tp{2S) radiative decays have been re- 
ported . BESII detector has accumulated 14 million 
i/>(2S) events, and it offers an opportunity to study a 
number of tp(2S) radiative decays. 

Experimentally, an event is reconstructed by us- 
ing the detector information on charged tracks and 
photons. For removing background events, the con- 
straints of the energy and momentum conservation 
for the final state of the signal channel are imposed 
on the fitting to the 4-component momentum (4C) 
reconstructed by the detector, and the background 
events can be suppressed by using information of the 
distribution of the \ 2 of the fit. With the help of 
the Monte-Carlo simulation, a reasonable requirement 
on x" 2 can be obtained by comparing the x 2 distribu- 
tions between the signal and the background chan- 
nels. However, the x 2 cut method fails badly some- 
times if the x 2 distribution of the background channel 
overlaps largely with the signal channels. In this pa- 
per, we try to fit the x 2 distribution to extract the 
number of signal events in measuring the exclusive 
decays of ip(2S) — > 7 + hadrons at BESII detector, 
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which includes 2, 4 and 6 prongs (stable hadron- 
pion,kaon, proton or anti-proton). 



II. FIT TO THE x 2 DISTRIBUTION 

Monte Carlo simulation shows that backgrounds of 
the radiative decay ofip(2S) — > 7+hadrons are largely 
coming from multi-photon decay channels, namely, 
the dominant background is ip(2S) — » 7r° + hadrons, 
and some contaminations are from ip(2S) — > 771"° + 
hadrons, along with other possible backgrounds. In 
principle, their x 2 distributions of the kinematic fit 
can be obtained with MCS. Since these types of back- 
grounds are due to missing photons, their distribu- 
tions of the x 2 shapes are different from each other. 
This property allows us to express the observed x 2 
distributions of the data sample by that of the signal 
channel and the background channels, i.e. 
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where w s and are the weight of the signal and the 
background decays, respectively. If the background 
channels are completely known, then the signal events 
can be extracted reliably by fitting to the data with 
Equation QJ. 



III. MONTE CARLO SIMULATION 

A. Background 

To make the measurement of ip' radiative decays re- 
liable, the background should be fully studied. With 
the help of SIMBES system [2j, extensive simulations 
have been performed. Table [I] summarizes the possi- 
ble background contributions to signal mode tp(2S) — * 
7 + hadrons. Generally speaking, the dominant back- 
grounds come from multi-photon decays, i.e. ip{2S) — » 
+ hadrons (n > 2). With more missing photons in 
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background decays, the contamination to the signal 
decay becomes lower and lower. 

Tab. I: The backgrounds for the radiative decay 
ip(2S) — ► 7 + hadrons, where the normalized back- 
ground events are based on 14M ip(2S) events and the 
branching fractions given in PDG table with require- 
ment m hadrons < 3.0 GeV. 



Signal mode 
i>(2S) - 


Background modes 


Nb g (normalized) 


1PP 


77T pp 


141.5 
8.8 


yK + K~n + n~ 




98 






16 


72(7r+7r _ ) 


7r u 2(7r+7r-) 


805.0 




TT + -K~ J/lfl, J/lj) — ► 3lT 


33.4 




■k"K + K-2(ty + tt-) 


6.0 




'yir°K + K-2(ir + iT-) 


6.1 



B. x 2 distribution 

Figure ^ shows the shapes of x 2 distributions for 
ip(2S) — > 7 + hadrons after data selection. In plots 
the points with error bars are the signal decays 
into 7pp (a), yK + K~tt + tt~ (6), 72(7r + 7r~) (c) and 
7i^ + i^~2(7r + 7r~) (d), respectively. The dashed his- 
tograms and the solid line histograms correspond to 
the distributions of the background channels ip{2S) — > 
7T° + hadrons and 777° + hadrons, respectively. For 
comparison of their shapes, the background events are 
respectively normalized to signal events. Obviously 
it can be seen that the shapes of the signal x 2 dis- 
tributions arc different from those of the background 
modes. This could be understandable due to the fact 
that the expectation of the x 2 value of the kinematic 
fit becomes larger to the backgrounds originated from 
one or more photons missing. 

C. Input-output check of \ 2 At method 

Since the dominant backgrounds of ip(2S) radiative 
decays originate from the photon missing decays, and 
their shapes are different from each other, therefore, 
fitting the x 2 distribution of the data with Equation 
provides us a tool to extract the signal events. 
We employ MCS to make an input-output check and 
estimate the uncertainty of the x 2 fit method. 

Figure [5] shows the x 2 distributions of the MC 
sample for ip(2S) decays into yK + K~TT + n~ (a), 
ir°K + K~tt + tt~ (b) and other unknown backgrounds 
by using 14M ip(2S) Lund-Charm MC sample. In 
order to make an input-output check of the x 2 fit 
method, a MC sample is generated with the x 2 shapes 
given in (a), (b) and (c) with an input number Nj? , 




Fig. 1: The chisqured distribution of MC 
samples for tp(2S) — > 7 + hadrons, where 
points with error bars show the \ 2 shapes for 
hadrons = pp (a), K + K~-k + tt~ (b), 2(-7r + 7r~) (c) 
and K + K~2(tt + -k~) (d), and the dashed histograms 
correspond to ip(2S) — ■» 7r° + hadrons and the solid line 
histograms to ip(2S) — * yir° + hadrons, respectively. 



iV^ and N£^, respectively. Then we extract the num- 
bers of events of these modes by the x 2 fit method. 
The total uncertainty includes errors from the x 2 fit 
method and the statistics, i.e. 



where is the number of the signal events obtained 
by the x 2 fit method, and cr st is the statistical uncer- 
tainty. Due to the uncertainty of the unknown back- 
ground shape, we also consider the unknown back- 
ground shape changed between the shape as shown in 
Figures [3 (b) and (c) parameterized by 

2 /-1 _ \ 2 , 2 
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where x 2 n an d X 2 g are the x 2 distribution of the 
unknown background and the dominant background, 
respectively. The x value is taken between zero to 
one. The MCS shows that if the x 2 shape of the 
unknown background tends to that of the dominant 
background, the uncertainty of extracting the number 
of the dominant background becomes large by x 2 fit 
method. However, this situation does not worsen the 
signal uncertainty. The MCS with a large statistics 
shows that the uncertainty of the x 2 fit method is less 
than 2%. 




x 2 
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IV. SUMMARY 

Based on MCS, it is found that the backgrounds of 
the ip(2S) decays into 7 + hadrons dominantly orig- 
inate from the missing photon decays like ip{2S) — > 
nj + hadrons (n > 2). The % 2 distribution shapes 
of backgrounds are distinctive from those of the sig- 
nal decays. This property allows us to extract the 
number of the signal events by the fit to the x 2 dis- 
tribution. The results of MCS study indicate that the 
uncertainty of the extracted signal events from the 
method is less than 2%. 



Fig. 2: The x 2 distribution for tp(2S) decays into 
the signal mode yK + K~ ir + ir~ (a), and the dominant 
background ty°K + K~it + iv~ (b) and other unknown 
background (c). In plot (d) the histogram corresponds 
to the x 2 distribution of the sample generated with 
X 2 shapes in (a), (b) and (c) and the solid line to the 
fitted results with these x? components. 
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